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ABSTRACT 

We have carried out a systematic search for the molecular ion C0+ in a sample of 8 
protoplanetary and planetary nebulae in order to determine the origin of the unexpect- 
edly strong HCO^ emission previously detected in these sources. An understanding of 
the HCO+ chemistry may provide direct clues to the physical and chemical evolution 
of planetary nebulae. We find that the integrated intensity of the C0+ line may be 
correlated with that of HCO+, suggesting that the reaction of CO^ with molecular 
hydrogen may be an important formation route for HCO+ in these planetary nebulae. 
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1 INTRODUCTION 

Planetary nebulae (PNe) are commonly identified and stud- 
ied by their optical emission arising in the ionized gas. 
However, CO and H2 surveys have revealed that a large 
fraction of th e gas in massive PNe is in fact in molecu- 



lar form (e.g. iHuggins et al.lll996l . l2005l: ISpeck et al 



2002 



1997 



Josselin fc Bachilleij |2003l . hereafter JB03) have shown that 



this molecular component may be dominant during a large 
part of the evolution of the nebulae and hence plays an im- 
portant role in the mass distribution and the shaping of the 
. PNe. The molecular gas in PNe will be mainly affected by 
■ the large radiation fields from the central stars: the dra- 
matic effects on the chemical composition will change as the 
star moves from the asymptotic giant branch (AGB) into 
' the protoplane tary nebula (PPN) phase and then into a PN 
(see review bv lKwok|[l993 ). 

Most chemical studies of the neutral gas in PNe have 
so far focused on comprehensive surveys of only a handful 
of PNe, all typically evolved, with the exception of CRL618 
and NGC7027 (e . g. ICox et al.1 Il992l: iBachiUer et all 1 19971 : 
iPardo et all 12004 120071 '). However. |JB03| recentlv reported 
high sensitivity observations of molecular lines in a sample 
of 7 young or intermediate-aged compact PNe, providing 
us with the best sample, so far, of PNe at different evolu- 
tionary stages. Their conclusions can be summarized as fol- 
lows: 1) some molecular species, such as HCO"'' and CN, are 
particularly abundant compared to envelopes around AGB 
stars or even to interstellar clouds; 2) a careful analysis of 
the chemical composition of a large enough sample can give 
an indication of the evolutionary stage of the PPN/PN; 3) 
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the survival of certain molecules can only occur in well pro- 
tected dense clumps (e.g. t he famous cometary knots in the 
Helix nebula and envelope, iHuggins et al.| [2002: Spec k et al.l 
[20021). The latter con c lusion supports the theoretical find- 
ings of iRedman et all l|2003l ). who showed that dumpiness 
in the slow wind from the central star may be responsible 
for t he high neut ral content of PNe, although other studies 
(e.g. IWoodsll2004] ) find that dumpiness in itself is not enough 
to protect the molecu lar c ompon ent. Interestingly however, 
iRedman et all l|2003l ) and |JB03| also agree in finding that 
HCO^ and HGN are severely underpredicted by theoretical 
models and, in fa ct, HCO^ is always fo und to be remarkably 
high in PNe (e.g. iBachiller et allflW?! ). 

Recent observ ational and theoretical work (e.g. JB03; 
I Woods et al.ll2003h has shown that the post- AGB phase is 
a period of great chemical changes. The appearance, or en- 
hancement, of HCO^ may correspond to different evolution- 
ary stages, depending on what causes its high abundance. 
For example, the standard path for HGO"*" formation, 

H;[ -l-CO ^ HCO+ +H2, (1) 

is generally excluded for PNe because it requires a high sur- 
vival rate of Hj[ , which is instead expected to be very short- 
lived due to its rapid dissociative recombination in the pres- 
ence of a UV radiation field. Conversely, harder radiation 
fields due to, for example, strongly ionizing X-ray emission 
from the central star work to maintain high abundances of 
HJ by increasing its production rate. If this is the cause of 
the high abundance of HCO"'' then the enhancement should 
occur at an early stage, just after the ionization of the nebula 
begins at the end of the PPN stage ('Kwok'l l993 ). An addi- 
tional formation route for HCO"'' is through the reaction of 
C0+ with Ha, 



CO-*- 4- Ha 



HCO+ -HH. 



(2) 
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Table 1. Properties of the observed sample of (proto-)planetary nebulae. 



Source 


KA(ZUUUj 
(hh mm ss) 


Dec. \Z\j\j\J) 

\ ) 


(kpcj 


bize 

(arcscc) 


V^LSR 
(km s "'" ) 


(km s "'" ) 


References 


BV5-1 


00 20 00.45 


+62 59 03.2 


5.5 


11 


-73 


10 


1,2 


K3-94 


03 36 08.09 


+60 03 46.3 


5.5 




-69 


16 


1,3 


CRL618 


04 42 53.67 


+36 06 53.2 


1.8 


14 


-22 


18 


4,5 


Ml-13 


07 21 14.95 


-18 08 36.9 


2.0 


<12 


+27 


18 


1,6 


Ml-17 


07 40 22.21 


-11 32 29.8 


2.5 


<12 


+28 


39 


1,7 


K3-24 


19 12 05.82 


+15 09 04.5 


3.5 


<io 


+44 


24 


1,3 


NGC 6781* 


19 18 28.09 


+06 32 19.3 


0.7 


110 


+17 


22 


4,8 


IC5117 


21 32 31.03 


+44 35 48.5 


3.0 


<12 


-11 


14 


1,6 



References: DlJosselin fc Bachil ieil ||2003|): 2')|josselin fc Bachillerl ll200lD: 3)^0886^ etHI 1I2OOOI'): 4) 



iBachiller et al.l lll997t)~5 ') Sanche z Contreras et al.l (|2004|) : 6) lHuggins et al.l l ll996l) : 7') lBachiller et al.l 
1 I1991I) : 8') iBachiller etal.1 819931 . 

* Since the angular size of this PN is significantly larger than the telescope beam, an offset position 
of (+10", +60") was adopted, at which HCO"*" had previously been observed. 



The reactive molecular ion CO can be formed by charge 
transfer between and CO or by direct ionization of CO 
by means of the propagation of a n ionizing front into molec- 
ular gas of high CO abundance (jLatter. Walker fc MalonevI 
1 19931 ). Thus, there are two formation routes for CO , one 
of which only becomes effective once the object is evolved 
enough for significant ionization of the nebula to take place 
(which occurs at late times). There are also two routes to 
form HCO"'", one that can only occur early in the evolution 
of the source and one that requires the presence of CO^. 
Therefore, the varying abundances of the two molecules can 
potentially be used as an evolutionary diagnostic tool. 

The obvious test to discriminate between the above sce- 
narios is to perform a systematic search for CO"'' in P Ne and 
correlate its emission with that of HCO , as, in fact, |JB03| 
suggest. CO"*" has b een clearly detected in one young PN, 
namely NGC 7027 dLatter et all Il993l : iHasegawa fc Kwokl 
I2OOII : iFuente et al.l |2003| ). In this paper we present the re- 
sults of such a search and show that it is likely that CO''' is 
directly correlated with HCO^, making these two ions good 
evolutionary tracers of PNe. 



2 OBSERVATIONS 
2.1 PNe sample 

We have selected a sample of 8 PPNe and PNe in which 
HCO^ has previously been detected. Our sample is primar- 
ily based on that of [jB0 3'. who detected the HCO'*' (1-0) 
line in 7 young to intermediate-age compact PNe. In ad- 
dition, we also include the PPN CRL 618 and the evolved 
PN NGC 6781 to increase the range of evolutionary stages 
in our sample. These two sources have also been pre- 
vious ly detected in HCO^(l-O) emission (IBachiller et al.l 
1 19971 ). We do n ot include NGC 7027, where CO+ has already 
been detected jLatter et al.l Il993l : iHaseeawa fc Kwokl l20oj : 
iFuente et al.ll2003l ). Table [Tl lists the properties of the PNe 
in our sample; their common names (column 1), J2000 coor- 
dinates of their central positions (2, 3), kinematic distances 
(4), based on the si mple approxim ation for the Galactic disc 
velocity structure of lBurtonI (|l974l ). sizes (5) based on CO(2- 
1) maps (where available; upper limits are given when the 
nebulae are unresolved), radial and expansion velocities, as 



measured in the C0(2-l) line (6, 7) and the papers from 
which this information was obtained (8). 



2.2 Observations 

The observations were obtained with the James Clarke 
Maxwell Telescope (JCMT), Mauna Kea, Hawaii on the 
nights of 2005 February 18-20 and 2005 October 2-10, with 
additional data collected on various dates in 2005. Het- 
erodyne receivers RxASi and RxB3 were used for the 200 
and 300 GHz band observations, respectively. Transitions of 
CO^ are described in terms of the quantum numbers A^, 
the rigid body angular momentum excluding electron spin, 
and J — N i: 1/2, the total angular momentum. The CO^ 
= 2 1, J = 5/2 ^ 3/2 (236.0626 GHz) line was 
searched for in the sample of 8 PNe, whilst detection of the 
CO'*' iV = 3 ^ 2, J = 7/2 ^ 5/2 (354.0142 GHz) line 
was only attempted in CRL 618, the source displaying the 
strongest emission in the lower CO^ line in our sample. The 
antenna half-power beamwidth (HPBW) was 20 arcsec at 
236 GHz and 14 arcsec at 354 GHz. A frequency resolution 
of 625 kHz was initially used with an instantaneous band- 
width of 920 MHz, however the resolution was changed to 
312.5 kHz and the bandwidth to 500 MHz for the data ob- 
tained in 2005 October. With the exception of NGC 6781, 
all observations were carried out in beam-switching mode, 
with a secondary mirror throw of 60 arcsec in azimuth . Since 
NGC 6781 has an angular size of 110 arcsec (Bachillc r et al.l 
Il993l ). the secondary mirror throw was increased to 180 arc- 
sec for this source to ensure that the nebula did not fall 
within the OFF position. Pointing and focus were checked 
every 60-90 minutes and the pointing uncertainty was found 
to be within 3 arcsec. The data were reduced in the stan- 
dard manner using the SPECX software package and were 
converted to the main beam (Tmb) scale using main beam ef- 
ficiencies determined by planetary observations made whilst 
the data were being acquired. During the data reduction pro- 
cess, several channels were binned and Hanning smoothing 
was applied to reduce the noise level. The channel binning 
was performed in the standard way, by averaging channels 
into bins 2-4 channels wide, starting from the first channel 
in the spectrum. 
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Figure 1. Spectra of the CO''" A' = 2 — > 1, J = 5/2 — > 3/2 line observed in the sample of PNe, with fits to the line profiles overlaid 
(dashed lines). Brightness temperatures are given in main beam (Tj^b) scale and velocities are quoted with respect to the local standard 
of rest (LSR). Some of the spectra have been binned to a lower velocity resolution to reduce the noise. 



3 RESULTS 

Spectra obtained for the sample of PNe are shown in Figs[T] 
and [21 We have fitted Gaussian profiles to the observed lines 
and these are overlaid on the spectra in Fig. [1] Table [2] lists 
the fit parameters for each source, including the integrated 
intensity (column 2), the peak main beam temperature (3), 
the RMS noise (4), the line width and central velocity (5, 
6), and the spectral resolution of the data (7, 8). For non- 



detections, upper limits for the line intensities are given by 

7(C0+) < 3a„,bAy, (3) 

where (Jmb is the main beam RMS noise of the spectrum and 
/Sy is the assumed line width, taken to be the expansion 
velocity of the source measured in the CO(2-l) line (see 
Table [T|. 

We obtain a strong (> IOct) detection of the 00"*" A'' — 
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Table 2. Summary of the Gaussian fit parameters for the observed CO"*" A'^ = 2 — » 1, J = 5/2 — > 3/2 line in each source. The integrated 
intensity, peak main beam temperature, RMS noise, line width, central velocity and the spectral resolution of the final data are listed 
for each source. In the case of non-detections, an upper limit to the integrated intensity is derived using equation |3] 



1 T*/^(21 


J ^ mb ^ ^ 






FWFTIVT 

1. VV i-±lVi- 


^SR 


Au 


Au 




(Kkms-1) 


(K) 


(K) 


(kms-1) 


(kms~^) 


(MHz) 


(kms-i) 


CRL618a 


0.762±0.077 


0.069 


0.008 


10.43 


-66.03 


1.25 


1.59 


CRL6186 


0.872±0.085 


0.062 


0.008 


13.24 


-24.81 


1.25 


1.59 


NGC 6781 


0.249±0.057 


0.031 


0.008 


7.60 


+8.73 


1.25 


1.59 


BV 5-1 


0.062±0.015 


0.011 


0.003 


5.18 


-63.82 


1.25 


1.59 


K3-24 


0.109±0.036 


0.013 


0.003 


8.16 


+39.86 


2.50 


3.17 


IC5117 


0.069±0.035 


0.010 


0.004 


6.22 


-10.64 


2.50 


3.17 


K3-94 


0.044±0.028 


0.010 


0.003 


4.17 


-72.41 


1.88 


2.38 


Ml-17 


<0.351 


0.006 


0.003 






2.50 


3.17 


Ml-13 


<0.162 


0.005 


0.003 






2.50 


3.17 



2 ^ 1, J = 5/2 ^ 3/2 line in CRL618 and marginal de- 
tections (^3(t) in another three PNe (NGC 6781, BV5-1, 
K 3-24). We obtain tentative detections in IC 5117 and K 3- 
94, and fail to detect the line in the remaining two PNe, 
M 1-13 and M 1-17, for which we derive upper limits to the 
line intensities. The linewidths of the previously detected 
HCO^ lines in these sources are not given in the literature, 
but are generally much wider than the 00"*" linewidths we 
find here. 

Estimates of the column densities have been derived 
under the assumptions of local thermodynamic equilibrium 
(LTE) and optically thin emi ssion, following the method 
outlined in lFuente et all (120031 ), 

iV = 1.94xlO«^^^gl)exp(^J^) I T^^dVxf^\ (4) 

where J Tmh dV is the integrated intensity of the CO^ line 
(Kkms^^ in main beam scale), /b is the beam filling fac- 
tor (described below), u is the frequency of the transition 
(in GHz), Q{Tcx) is the partition function for CO^ at the 
excitation temperature Tex, evaluated by interpolating be- 
tween values obtained from the Jet Propulsion Lab oratory 
(JPL) molecular line data base jPickett et al.|[l998l ). Eu is 
the upper state energy, is the upper state degeneracy and 
Aij is the Einstein spontaneous emission coefficient for the 
transition (in s~^). The relevant values for the A*' = 2 ^ 1, 
J = 5/2 — + 3/2 line are listed in Table |3] The beam filling 
factor is given by 

where 6's is the source size (FHWM) and 6h is the telescope 
beamwidth (HPBW) . For the source sizes, we have adopted 
the CO(2-l) measurements listed in Table [l] We note that 
the 00"*" emission may arise from smaller regions of the 
nebulae than are traced by CO, however, there are no maps 
of CO'*' (and very few of HCO^) that may be used, whereas 
maps of CO emission exist for most of the sources. When 
sources have not been resolved in previous observations, we 
have used the upper limits on their size. 

Since we would like to correlate the CO"*" emission with 
that of HCO+, we have assumed that the two species are 
emitted from the same region and we have adopted the 
sa me uni form excitation temperature Tox = 25 K assumed 
by |JB03| in their calculation of the column densities of 
HCO^, with the exception of CRL618, for which we de- 
tected two 00"*" transitions and can therefore derive an es- 



timate for the excitation temperature using the rotation di- 
agram method (see Sect. 13. ip . The chosen temperature of 
25 K has been deduced from CO observations in several neb- 
ulae l|Bachiller et all 119931 ). As noted before, however, CO 
is probably ubiquitously present in the molecular gas sur- 
rounding these objects, whereas HCO"*" and CO^ may only 
exist in the thin layers of photodissociation regions (PDRs), 
hence we may be underestimating the kinetic temperature. 
When the line was not detected, upper limits determined 
for the integrated intensities were used to derive upper lim- 
its for the column densities. The column densities are re- 
ported in Table [3] where, for reference, we also list those 
derived previously for HCO"'', again corrected for the filling 
factor (assuming the same source sizes and an TRAM 30-m 
beamwidth of 28 arcsec). It should be noted that the uncer- 
tainties quoted in the table are likely to be overshadowed 
by larger systematic errors, including uncertainties in the 
beam filling factors, and may therefore be misleading when 
comparing the column densities. 



3.1 CRL618 

The PPN CRL 618 displays the strongest emission in our 
sample and has been clearly detected in two lines of CO^ 
{N = 2 ~* I, J = 5/2 ^ 3/2 and iV = 3 ^ 2, 
J = 7/2 5/2; see Fig. [2|. This multiple line detection al- 
lows the rotational temperature of the gas to be determined 
with a rota tional (Boltzmann) diagram using the method 
detailed by iFuente et all (|2003l ). A straight line fit to the 
rotational diagram yields a rotational temperature of 32 K. 
This higher temperature is si milar to the assu med excita- 
tion temperature of 25 K CBac hiller et al.l | l993l ), but lower 
than would be expected for a PDR region. If taken to be 
representative of the temperature for all of our sample, an 
excitation temperature of 32 K reduces the column densities 
presented in Table |4] by 14 per cent. 

Two velocity components are clearly seen at —70 and 
—30 kms~^ in Fig. [51 in both transitions of CO^. These ve- 
locities are consistent w ith the interferometric observation s 
of HCO+ and HCN bv lSanchez Contreras fc Sahail l|2004h . 
They detected HCO+ J = 1 ^ and HCN J = 1 ^ emis- 
sion across a velocity range —120 Vlsr ^ +40 kms~^, 
arising from two distinct regions, which they assigned to a 
high-velocity bipolar outflow and a slowly expanding equa- 
torial torus. We identify the two blue-shifted features in 
the CO"'' transitions as follows. Feature a is the fast com- 
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Figure 2. Spectra of the CO+ Af = 2 1, J = 5/2 -» 3/2 
and AT = 3 ^ 2, J = 7/2 ^ 5/2 lines observed in CRL618. 
Brightness temperatures are given in main beam (Tj^b) scale and 
velocities are quoted with respect to the local standard of rest 
(LSR). The two velocity components appearing in each of the 
lines are labelled a and b on the spectra. 



ponent, arising in the blue-shifted (i.e. brightest) lobe of 
the bipolar outflow, and feature 6 is due to the slowly ex- 
panding torus. We propose that the red-shifted lobe of the 
bipolar outflow is not visible in CO^ emission because it 
is much f ainter, as is the case for the HCQ" *" emission ob- 
served by ISanchez Contreras fc Sahail l|2004l ) . Their fluxes 
were 2.7 (red) versus 6.4 (blue) Jy beam~^ kms~^. The slight 
offset in peak velocities in the two transitions in Fig. [2] is 
likely due to differing e xcitation conditions in the gas; the 
HCO+ and HON lines in lSanchez Contreras fc Sahail (|2004l ) 
also peak at slightly different velocities. Interestingly, our 
CO^ line components are of comparable strengths. This 
contrasts with the HCO'*' and HCN lines, in which the 
outflo w component is much weaker tha n the torus compo- 
nent jSanchez Contreras fc Sahail |2003 . their flgure 1), and 
may indicate that different CO"*" formation and destruction 
routes operate in these two regions. 



4 HCO+ AS AN EVOLUTIONARY TRACER? 

It is impossible at this stage to establish if and when CO^ is 
the source of the HCO^ enhancement. Nevertheless, despite 
the lack of a statistically meaningful sample, in Figs[3]and|3] 
we plot the derived column densities of CO^ and HCO^ (the 
latter taken from the literature) for our sample of planetary 



Parameter 


AT = 2 ^ 1, 
J = 5/2 ^ 3/2 


Af = 3 ^ 2, 
J = 7/2-* 5/2 


y (GHz) 


236.0626 


354.0142 


Eu (K) 


16.997 


33.987 


/I (Debye) 


2.771 


2.771 


Q(25K) 


18.35 


18.35 


gu 


6 


8 


A,j (s-1) 


4.7x10"* 


1.7x10-3 



Table 4. Estimated column densities of CO^ 
and HCO+. 



Source 



Af(CO+) 
(cm-^) 



CRL 618 

NGC 7027 

Ml-17 

IC5117 

Ml-13 

BV5-1 

K3-94 

K3-24 

NGC 6781 



3.7±0.4 
1.8±0.4 

<1 
3.6±1 

<8.5 
3.7±0.9 
3.1±2.0 
7.6±2.5 
3.5±0.8 



xlO^^ 
xlO^^t 

X1012 
xlO" 
xlO" 
xlO" 
xlO" 
xlO" 
xlO" 



Af(HCO+) 
(cm-^) 
IxlO^-** 
lxlO"t 
IxlOiat 
6x10^21 
6xl0i2t 
2x10^21 
3x10^21 
3xl0i2t 
5x10^2* 



1.6±0 
4.3±0 
l.litO 
3.5±0 
9.0±0 
1.3±0 
2.3±0 
5.3±0 
5.5±0 



* Values taken from 
t Values taken from 



Bachiller et al.l lll997l) 



Josselin fc BachilleJ 1 120031) 



Calculated using the CO^ int egrated intensity 
observed bv lLatter et al.l lll993l) 



nebulae, i ncluding NGC 702 7 (for which we take the CO 
data from lLatter et aOi 19931 and HCO"'' data from jjBO^ ) as 
a function of their believed evolutionary stage, from PPNe 
(left) to evolved PNe (right), as well as their ratios. The 
evolutionary stage of each PNe is estimated based on the 
ratio of the molecular mass to the ionized mass Mm/Mi, 
which is believed to be a good indicator and almost inde- 
pendent of distance (iHuggi ns et al.l[l996[ l. Fig. [3] does seem 
to indicate a correlation between HCO'*' and CO^. For ex- 
ample, both CRL 618 and NGC 7027, believed to be young 
objects, have the highest HCO^ and the highest CO^ abun- 
dances (in fact they are the only 2 objects for which a strong 
detection of CO^ has been reported). This implies that, al- 
though young and hence enough should still be available 
to form HCO^, CO^ is indeed available and, at early evolu- 
tionary stages, must be formed by charge transfer between 
H"*" and CO. It is then also interesting to note that as we 
move along the evolutionary a;-axis, after a dip, both HCO"'' 
and CO"'' tend to increase with age. This may mean that at 
later stages, when the propagation front appears, more CO^ 
(and hence 1100"*") is formed (via the second route of forma- 
tion). We cannot draw any conclusions from the 2 objects 
for which CO^ was not detected (for which we derive an 
upper limit) but for which HCO^ is quite high. A possible 
explanation for the lack of CO^ is that at this evolutionary 
stage the radiation is strong enough to destroy this species, 
whilst there is still enough to make HCO'*'. 



5 CONCLUSIONS 

We have searched for the molecular ion CO^ in a sample 
of PPNe and PNe where HCO"'' was previously detected. 
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Figure 3. Column densities (with error-bars) of CO"*" (diamonds) 
and HCO"*" (crosses, from tiie literature) for the sample of plan- 
etary nebulae. The sources span a range of evolutionary stages 
from PPN (left) through to evolved PN (right). Only an upper 
limit was obtained for the column density of CO"*" for M 1-17 and 
M 1-13. 
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Figure 4. Ratios of CO"*" to HCO^ column density (with error- 
bars) for the sample of planetary nebulae. The sources span a 
range of evolutionary stages from PPN (left) through to evolved 
PN (right). Only an upper limit was obtained for the column 
density of CO+ for M 1-17 and M 1-13. 



If detected, CO^ may help our understanding of the origin 
of the strong HCO^ emission. We obtain a strong detec- 
tion of CO^ in two transitions in CRL618 and > 3cr detec- 
tions of the lower transition in three other objects, namely 
NGC 6781, BV5-1 and K3-24. The lower CO+ transition 
is also weakly detected in IC5117 and K3-94, but is not 
detected in M 1-13 and M 1-17. We find that the integrated 
intensity of the 00"*" line is reasonably well correlated with 
that of HCO^. Although we cannot draw any strong con- 
clusion on the origin of the HCO"'" emission, we can con- 
clude that CO^ is present in at least some early PNe and 
hence this may provide a supplementary source for HCO"'' 
and that, in evolved PNe, where Hjj" should be almost ex- 
hausted, 00"*" may be present in high enough abundances 
to boost HCO^ again. 
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